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Abstract 

We compute the running QCD coupling on the lattice by evalu- 
ating two-point and three-point off-shell gluon Green's functions in 
a fixed gauge and imposing non-perturbative renormalisation condi- 
tions on them. Our exploratory study is performed in the quenched 
approximation at j3 = 6.0 on 16 4 and 24 4 lattices. We show that, 
for momenta in the range 1.8 — 2.3 GeV, our coupling runs according 
to the two-loop asymptotic formula, allowing a precise determination 
of the corresponding A parameter. The role of lattice artifacts and 
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Studi di Milano, Via Celoria 16, 20133-Milano, Italy. 
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finite-volume effects is carefully analysed and these appear to be un- 
der control in the momentum range of interest. Our renormalisation 
procedure corresponds to a momentum subtraction scheme in con- 
tinuum field theory, and therefore lattice perturbation theory is not 
needed in order to match our results to the MS scheme, thus elim- 
inating a major source of uncertainty in the determination of a^. 
Our method can be applied directly to the unquenched case. 
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1 Introduction 



The running coupling a s (/i), where fi is a momentum scale, is a fundamental 
QCD quantity providing the link between low and high-energy properties of 
the theory. Given a renormalisation scheme, a s (fi) can be measured exper- 
imentally for a wide range of momenta. A precise determination of a s (/x) 
(or equivalently of the scale A determining the rate at which a s runs) is ex- 
tremely important as it would fix the value of a fundamental parameter in 
the Standard Model, providing bounds on new physics. 

Computing a s is a major challenge for the lattice community. Several dif- 
ferent lattice definitions of the renormalised coupling have so far been inves- 
tigated [[[], |2], ||, f|. Apart from the static quark potential approach ||, which 
involves a phenomenological parametrisation of the interquark potential, one 
feature of all other definitions of the coupling is that lattice perturbation the- 
ory (LPTH) has been used in order to convert the measured numerical value 
into a value for cttfs- Despite recent proposals to improve the convergence of 
LPTH series ||, this step still provides one important source of systematic 
errors in the final prediction for a-Ms- At present, systematic errors, namely 
quenching, discretisation effects, finite volume effects and LPTH, dominate 
statistical ones for state-of-the-art computations. 

We investigate here a more recent proposal |J for the determination of 
a s from the renormalised three-gluon vertex. This is achieved by evaluating 
two-point and three-point off-shell Green's functions of the gluon field on the 
lattice, in the Landau gauge, and imposing non-perturbative renormalisation 
conditions on them, for different values of the external momenta. By varying 
the renormalisation scale /i, one can determine a s (ii) for different momenta 
from a single simulation and analyse the /^-dependence of the coupling. In 
particular, one can investigate if the asymptotic behaviour is reproduced for 
large momenta. In practice, for a given choice of the lattice parameters, one 
needs to choose fi in a range of lattice momenta such that both finite volume 
effects and discretisation errors are under control. If such a momentum 
region exists and if the coupling is found to run according to the two-loop 
asymptotic formula, then we can get a meaningful measurement of 
in our renormalisation scheme which can then be related perturbatively to 
other definitions of the coupling. 
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As will be explained in the following, a crucial feature of this procedure, 
which corresponds to momentum subtraction renormalisation in continuum 
QCD, is that renormalised Green's functions do not depend on the way the 
theory is regularised. As a consequence, LPTH is not needed in order to relate 
the measured coupling to a^, and the relation between the two schemes can 
be computed entirely in continuum perturbation theory. An analogous non- 
perturbative approach has been recently applied, with encouraging results, 
to the renormalisation of composite fermion operators ||. 

The aim of this paper is to demonstrate the feasibility of our approach 
and to investigate the role of systematic lattice uncertainties, such as dis- 
cretisation effects and volume dependence. We perform this investigation 
in the quenched approximation at (3 — 6.0 for two different lattice volumes. 
Given the simplicity of the method, we expect the application to full QCD 
to present no additional problems 0. 

The paper is organised as follows: in Section |2| we present the renormal- 
isation scheme, the definition of the coupling and the numerical procedure. 
In Section |3] we give the lattice results, discussing all sources of systematic 
errors in the calculation. In Section |] we discuss the procedure to relate our 
coupling to ajfs and we compare our non-perturbative results with predic- 
tions based on the use of LPTH. In Section [j| we summarise our work and 
we discuss some possible future developments. Finally, the details of the 
perturbative calculations are given in the appendix. 



2 The Method 

Provided systematic errors are under control, in order to compute a-^s from 
the lattice it is sufficient to: 



1. set the scale of momenta in physical units by determining the lattice 
spacing a; 

2. define a suitable renormalisation scheme and a renormalised coupling 
to be measured; 



3. match the result to the MS scheme. 
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The physical quantities most frequently used to determine the value of a 
are: the p meson mass, the string tension, the IP-IS mass splitting in heavy 
quarkonia |]J, and a characteristic length r , phenomenologically connected 
to the intermediate range of the heavy quark potential ||. Each choice has its 
theoretical and technical advantages, extensively discussed in the literature. 

In this work we set (3 = 6.0 and we take the value of a" 1 determined 
by Bali and Schilling |§ in their string tension measurements. These yield 
a -1 = 1.9 ± 0.1 GeV. We quote a systematic error on the scale to take into 
account the uncertainty resulting from other possible choices. 

The definition and measurement of the coupling is achieved by computing 
the gluon propagator and the three-gluon vertex function in a fixed gauge 
and renormalising them in a non-perturbative way. We define the lattice 
gluon field A^(x) as 



2iag 

where fi indicates a unit lattice vector in the fi direction and go is the bare 
coupling constant (we omit the colour index). 

After performing the Fourier transform of ([!]) one can define (unrenor- 
malised) lattice n-point gluon Green's functions, in momentum space: 

Gu\ W2 ...^(pi,P2,---,Pn) = (A^ip^A^^.-.A^ipn)), (2) 

where (■) indicates the Monte-Carlo average and momentum conservation 
implies p\ + p 2 + ■ ■ ■ + p n = 0. 

Since the lattice calculation aims to evaluate such Green's functions in 
the "continuum window", i.e. for a range of parameters such that continuum 
physics is observed, the following criteria must be satisfied: 

1. (3 > 6.0, so that scaling is observed for physical quantities. 

2. La is large enough in physical units, where L is the linear dimension 
of the lattice, such that finite-volume effects are under control. 

3. Discretisation errors due to the contribution of nonleading terms in a 2 
to the Green's functions are negligible. 



5 



Since we work in momentum space, the last requirement can be directly 
translated into a restriction on the range of lattice momenta that should 
be used. In the remainder of this section we assume that all the above 
requirements are met and we adopt the formalism of continuum QCD. 

We choose to work in the Landau gauge. As already mentioned, the 
quantities of interest are the unrenormalised gluon propagator 

G [ u\Ap) = T» v (p) Gu(p 2 ), TM ^hv- Vj f- (3) 

.(3) 



and the complete unrenormalised gluon three-point function G v Q/ g 7 (pi, P2, Pz)- 



This is related to the one-particle- irreducible three-point function Tjj a/3 (pi,p 2 , 



by the formula 

G !/V(Pl)P2,P3) = r { U Sp z(pi,P2,P3) G^Sa.ipi) GP pl3 (j>2) G^ 7 (p 3 ). (4) 



In order to renormalise the gluon wave function, we impose that at a 
fixed momentum scale p 2 = fi 2 the renormalised gluon propagator takes its 
continuum tree-level value. One gets the non-perturbative renormalisation 
condition 

G R (p)\ p2= ^ = Z2 1 (jMi)G u (j)a)\ t fl =l jfi = — . (5) 

/i 

The above equation defines in a non-perturbative way the gluon wave-function 
renormalisation Za- 

For the three-gluon vertex, we start by choosing a kinematics suitable for 
the lattice geometry and which allows a simple non-perturbative definition 

(3) 

of the vertex renormalisation constant Zy. Recalling the general form of Y\j 
in the continuum |10[, it turns out that if one evaluates G v in the Landau 



gauge at the kinematical points defined by 

a = 7» Pi = -Pa =P, P2 = 0, (6) 

then one can write 



E«=i Q/ Qj g a (pa,0,-pa) 



6 i Z v 1 (pa) g p p . (7) 



(Gu(pa)y Gu(0) 

With the above definition Zy contains a term which is linear in the external 
momenta but not proportional to the tree-level vertex (see appendix). Notice 
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that the numerical calculation on the lattice for the three-gluon vertex yields 
directly the product Zy^g^. 

At this point we can define the running coupling g at the scale \i from 
the renormalised three-gluon vertex at the asymmetric point as 

9(p) = z Tin a ) Zy X {jm) g , (8) 

where the relevant renormalisation constants have been defined in (|5]), (|7|) 
and a s (n) = g(/j,) 2 /A7c. This choice corresponds to a momentum subtraction 
scheme, usually referred to as MOM in continuum QCD fllf . We postpone 
the discussion of the matching procedure until Section f|. 



2.1 Computational Procedure 

SU (3) gauge configurations were generated at (3 = 6.0 at two lattice sizes; 150 
configurations on a 16 4 lattice, and 103 configurations on a 24 4 lattice. The 
configurations on the smaller lattice were generated on a 16K CM-200 at the 
University of Edinburgh, using a hybrid-overrelaxed algorithm, where both 
Cabibbo-Marinari pseudo-heatbath and overrelaxed updates were performed 
on three SU(2) subgroups. Successive configurations were separated by 150 
sweeps (every sixth sweep a heat-bath), 1000 sweeps being allowed for ther- 
malisation. Landau gauge-fixing was achieved by using a Fourier-accelerated 
algorithm . Autocorrelations were investigated by performing a standard 
jackknife error analysis. 

The data for the larger lattice were generated on the ACPMAPS super- 
computer at FNAL using the Creutz pseudo-heatbath algorithm, with 1600 
sweeps between configurations. The configurations were fixed to the Lan- 
dau gauge using an overrelaxation algorithm, with the final iterations being 
performed in double precision. 

A crucial step in the method is the accurate implementation of the lattice 
Landau gauge condition 

A(z) =J2M x + ft-M x ) = Q - ( 9 ) 

To monitor the gauge-fixing accuracy we compute the quantity 

d = TArT,^ At(x)A(x) (10) 

v C x 
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as the algorithm progresses, terminating when 9 < 10~ n (this is close to 32- 
bit machine precision). Since our calculations involve low momentum modes 
of gluon correlation functions, the above test is not in itself sufficient because 
A(x) is a local quantity. For this reason, we also compute 

A (t)=X>o(:?,t). (11) 

x 

In a periodic box, the Landau gauge condition implies that A (t) is inde- 



pendent of t JT3|]. In particular, for our configurations on the smaller lattice, 
A (t) was constant to better than one part in 10 5 . For the purpose of our 
analysis, the only quantity that needs to be stored is the Fourier-transformed 
field A^{jp) for a selected range of lattice momenta. All n-point gluon cor- 
relation functions can then be assembled using eq. @, where momentum 
conservation is imposed explicitly. 

To compute a s , we first evaluate the gluon propagator and determine Za 
from eq. (||]). Next, we measure the complete three-point function G$ of the 
gluon field and the quantity on the l.h.s. of (|7|). Finally, g(/j,) is obtained from 
eq. (§). We take advantage of all the symmetries of the problem to improve 
statistics. The quoted errors are obtained using a single-point-elimination 
jackknife algorithm. 



3 Results 

3.1 Tensor Structure 

We start by analysing the tensor structure of the lattice gluon propagator 
and three-gluon vertex function as a means both of determining the degree 
of violation of continuum rotational invariance, and of verifying the extent 
to which the Landau gauge condition is satisfied in momentum space. 

It is worth noting that our definition QXj) for the gluon field differs from 
the one which has been used in all non-perturbative calculations to date, 
which is ||13||: 



t/,.w-^w_i Ti /^)-^wy 



2iago 3 V 2iag, 
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It turns out that the above "asymmetric" definition is not consistent with 
the one which is usually used in perturbative lattice calculations. Consis- 
tency is achieved by using the "symmetric" definition ([I]), which has better 
properties in the continuum limit. To illustrate this point, we observe that 
in momentum space the two definitions are related by the formula 

AM = e-^ /2 A». (13) 

If we now write the Landau gauge-fixing condition (^) in momentum space, 
using the asymmetric definition one gets 

53 [(cos^v - 1) - i sinp M ] A^(p) = 0, (14) 

while the symmetric definition yields 

£2isinp M /2A» = 0. (15) 

In the limit a — ► the continuum gauge condition is recovered with 0(a) 
corrections in the asymmetric case and 0(a 2 ) in the symmetric one. Thus 
the latter corresponds to an "improved" lattice Landau gauge condition. 

By using the symmetric definition one can check very accurately the 
tensor structure of non-perturbative lattice Green's functions against what 
is expected from LPTH. Based on such a definition we expect the Landau 
gauge propagator to satisfy 

G { u\M = %M Gu(p 2 ), f,M = V - ^ (16) 

where 

This makes the analysis of violations of rotational invariance very simple 
since, if our propagator is found to satisfy (0), one has for any p ^ 0: 

E Gf^p) = 3 Gu(p% (18) 

so that violations of rotational invariance in G v (p) can only arise 
from the scalar part Gu(p 2 ). If the asymmetric definition is used, the simple 
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V 


Po sinpo/2 Goi 
pi sinpi/2 Goo 


pi sinpi/2 Goi 
po sinpo/2 Gn 


p'f sin 2 pi/2 Goo 
pj, sin^po/2 Gn 


(1,1,0,0) 
(1,2,0,0) 
(2,2,0,0) 
(1,3,0,0) 
(3,3,0,0) 


1 1 1.00 
1/2 0.510 0.510 

1 1 1.00 
1/3 0.351 0.351 

1 1 1.00 


1 1.0 1.00 

2 1.962 1.962 
1 1.0 1.00 

3 2.848 2.848 
1 1.0 1.00 


1 1.0 1.00 
4 3.849 3.848 
1 1.0 1.00 
9 8.110 8.110 
1 1.0 1.00 



Table 1: Symmetry tests for G v on the 16 4 lattices at j3 = 6.0, using 
the symmetric definition of A^. 



p 


po sinpo/2 Goi 
Pi sinpi/2 Goo 


Pl sinpi/2 Goi 
Po sinpo/2 Gn 


p{ sin 2 pi/2 Goo 
p^ sin 2 p /2 Gn 


(1,1,0,0) 
(1,2,0,0) 
(2,2,0,0) 
(1,3,0,0) 
(3,3,0,0) 


1 1 1.00 
1/2 0.510 0.50 

1 1 1.00 
1/3 0.351 0.324 

1 1 1.00 


1 1 1.00 

2 1.962 1.924 
1 1 1.00 

3 2.848 2.631 
1 1 1.00 


1 1 1.00 
4 3.848 3.848 
1 1 1.00 
9 8.110 8.110 
1 1 1.00 



Table 2: Symmetry tests for Gy (jp) on the 16 4 lattices at (3 — 6.0, using 
the asymmetric definition of A^. 



tensor structure in eq. (16) is modified and the interpretation of the numerical 
results is more difficult. 

In Tables [I] and we show ratios of tensor components of the gluon prop- 
agator on the 16 4 lattices using the symmetric and asymmetric definitions of 
A^ respectively. These are compared to what is expected from ( |TED and its 
continuum counterpart. In all cases the uncertainty is less than one unit in 
the last quoted figure. We note the following: 

• For the symmetric definition of A^, the numerical data are completely 
consistent with our expectation from eq. (|16D . This indicates that the 
lattice gauge condition has been implemented very accurately, so that 
we have a complete understanding of the tensor structure of the two- 
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point function. 



• For the asymmetric definition of A^, the picture is quite different, al- 
though even in this case the deviation from (|T6|) is at most 10% for the 
momenta shown. 

We consider now the three- gluon correlator G v a ^{p, 0, —p). On the lat- 
tice we expect it to satisfy 

Gf opyfo 0, -p) = f aj (p) F(p 2 ) p p (Gu(p 2 )) 2 Gu(0), (19) 

which is the form one obtains by performing the substitution p^ — > p^ in 
the continuum expression. Here F(p 2 ) represents a generic function of p 2 . 
Note that from (|l3j ) it follows that Gu a g (p, 0, —p) does not depend on the 
definition used for the gauge field as the extra phase factor associated with 
the symmetric definition of A^(p) cancels at this kinematic point. 

In Tables | and f| we show ratios of tensor components of G^ a g (p, 0, —p) 
on the 16 4 and 24 4 lattices compared to expectations from (|19D and its con- 
tinuum counterpart. 



p 


p{ sin 2 pi/2 Goio 
pj, sin 2 p /2 Gm 


PO sinpo/2 Gioi 
pi sinpi/2 Gm 


pi sinpi/2 Gon 
po sinpo/2 Gm 


(1,1,0,0) 
(1,2,0,0) 
(2,2,0,0) 


1 1 1.000 
4 3.848 3.848 
1 1 1.000 


1 1 1.1(2) 
1/2 0.510 0.3(1) 
1 1 0.8(3) 


1 1 1.000 

2 1.962 1.962 
1 1 1.000 



Table 3: Symmetry tests for G v a/ g 7 (p, 0, —p) in the Landau gauge, using the 
symmetric definition of A^, on 16 4 lattices. Unless otherwise noted, the error 
is always less than one unit in the last quoted figure. 

The theoretical expectation is satisfied very accurately for two of the 
three ratios under consideration, while the agreement for the ratio Gaoi/Cm 
is much poorer, especially on the larger lattice. This is related to the fact that 
the first two ratios are completely determined by the Landau gauge condition 
for our choice of the momenta, whilst the last ratio is not. Thus these 
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p 


p( sin 2 pi/2 Goio 


Po sinpo/2 Gioi 


P! sinpi/2 Gnu 


sin 2 p /2 Gm 


pi sinpi/2 Gm 


po sinpo/2 Gm 


(1,1,0,0) 
(1,2,0,0) 
(2,2,0,0) 
(1,3,0,0) 


1 1 1.000 
4 3.932 3.932 
1 1 1.000 
9 8.596 8.596 


1 1 -0.3(4) 
1/2 0.504 -0.2(4) 

1 1 0.3(7) 
1/3 0.341 0.8(7) 


1 1 1.000 

2 1.983 1.983 
1 1 1.000 

3 2.932 2.932 



Table 4: Symmetry tests for G\j a g 7 (p, 0, — p) in the Landau gauge, using the 
symmetric definition of A^, on 24 4 lattices. Unless otherwise noted, the error 
is always less than one unit in the last quoted figure. 



results provide further evidence that the momentum-space Landau gauge- 
fixing condition is very well satisfied on all our lattices, but they also hint at 
a poorer quality of the data on the larger lattice. 

3.2 Renormalisation Constants and Running Coupling 

We are now left with the task of computing the renormalisation constants 
defined in (|5|) and (|7|) and the running coupling (||). Notice that because of 
our definitions and choice of kinematics, the difference between the symmetric 
and asymmetric definitions for the gauge field is immaterial for this purpose. 
All data are plotted vs. /i = y/p 1 , expressed in GeV. In order to detect 
violations of rotational invariance, we have used whenever possible different 
combinations of lattice vectors for a fixed value of p 2 and we have plotted 
separately the corresponding data points. 

In Figure [I] we show the scalar gluon self-energy Gu(fi). Statistical errors 
are negligible and no violations of rotational invariance can be detected in the 
momentum range under consideration. Next we compute the gluon three- 
point function and evaluate the running coupling according to (H). This is 
plotted in Figure |2| for both lattice sizes. We obtain again a clear signal; on 
the other hand, the data do show some violations of rotational invariance. 
We note again that the errors on the 24 4 lattice are considerably larger than 
on the smaller lattice. 
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Figure 1: vs. ji for the 16 4 lattices at (3 = 6.0. 

3.3 Systematic Lattice Uncertainties and Extraction 
of A 

The important physical question is whether one can isolate a range of mo- 
menta where lattice artifacts are negligible and our coupling runs according 
to the two-loop perturbative expression 



9*M 



(20) 



where b = 11/167T 2 , b\ = 102/(167r 2 ) 2 and J^^y M is the QCD scale pa- 
rameter for the renormalisation scheme that we are using (in the quenched 
approximation). To answer this question, and obtain an estimate for A^^g^., 
we compute A JV ~ 7 . as a function of the measured values of g 2 (/i) according 
to the formula 

A MOM = ^p(-^))Mrt]"^. (21) 

If the coupling runs according to ( p0|) , then ^ A ^f M as defined from the above 
equation must be constant. Given the exponential dependence of ^j^y M on 
g 2 (/i), this test is a very stringent one. As an alternative procedure one could 
fit the data for g 2 (n) to formula (|20"D , with ^^q M as a parameter, but this 
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Figure 2: Running coupling vs. /i. 
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1.5 



1.0 — 



< 
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0.0 



Figure 3: ^j^f M vs. fi for our two lattice sizes. 



would not result in a strong test since (|20|) only depends logarithmically on 

MOM' 

By plotting A^^^. versus \i (see Figure 3), three different regimes can 
be identified: 



1. For fM < 1.8 GeV, ^ h ^=f M displays a strong dependence on the renor- 
malisation scale. This is not surprising, as for low momenta asymptotic 
scaling is not expected. 

2. In the range 1.8 < /i < 2.3 GeV the data are consistent with a constant 
value for A A ^^ f . No violations of rotational invariance are observed in 
such a range and a comparison of the two lattice sizes shows no volume 
dependence either. This is shown in Figure (|. 

3. For fi > 2.3 GeV, rotational invariance is broken by higher order terms 
in a 2 and the two-loop behaviour disappears. 



In summary, we appear to have a "continuum window" in the range 
1.8 < fi < 2.3 GeV, where two-loop scaling is observed and lattice artifacts 
are under control. In order to extract a prediction for A^- M , we fit the data 
points in the continuum window to a constant. We take as our best estimate 
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1.4 



1.2 



16 



< 



0.8 



24 



1.0 



2.0 2.2 



Figure 4: ^^^ M vs. /i in the "continuum window" . 
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the fit to the 16 4 data, for which the statistical errors are smaller, and obtain 



A~ m = 0.96 ± 0.02 ± 0.09 GeV, (22) 

where the first error is statistical and the second error comes from the un- 
certainty on the value of a" 1 . 

Since our lattice calculation is performed in a fixed gauge, it is worth 
mentioning that non-perturbative gauge-fixing ambiguities (Gribov copies) 
may in principle affect our results. These could be regarded as a potential 
source of systematic errors. However, based on previous investigations of 



other lattice gauge-fixed quantities[14], we believe that the contribution of 



'Gribov noise" to the overall error, if at all present, is negligible. 



4 Matching to MS 



In this section we extract a prediction for ocjjg with zero active quark flavours 
from our numerical results for ^j^ M - As already mentioned, the procedure 
that we adopt avoids the use of LPTH. We also discuss a LPTH calculation 
which provides some important consistency checks. 

We start from our numerical estimate (p2|) of the continuum scale pa- 
rameter h-jtfQ M - The general relation between the scale parameters in two 
continuum schemes A and B can be written as 



A 



A 

A^ = eXP 



1 / 1 



with 



)+0(g 2 ^))}, (23) 



M = 9bM{i + ^Cab + 0(g 2 M), (24) 

where Cab, which in general depends on the chosen gauge and on the number 
of active flavours, is obtained from a continuum perturbative calculation as 
described in the appendix. Since A^, A# are independent of fi and because 
of asymptotic freedom, by letting fi —>■ oo the ratio (|23"D can be determined 
to all orders in the coupling constant from the one- loop calculation (see for 
example [15|, [T6|j). 
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We obtain, in the Landau gauge and for zero quark flavours 



^ = 0.35. (25) 

MOM 

This result is in agreement with with previous one-loop calculations of the 
three-gluon vertex [I7|. 



Note that the scheme that we have adopted for the non-perturbative 
calculation (see Section 2) differs from the usual MOM scheme as it contains 
an extra constant term in the vertex renormalisation constant. This is a 
linear term in the momenta, not proportional to the tree level vertex, which 
is equal on the lattice and in the continuum. The coefficient in ([25]) includes 
it perturbatively (see appendix). 

Using (p2] ) and (|25|) , we can extract from 



MS 
a(0) 

Ag> =A ( °^ -^F-- (26) 

A/5 MOM a(0) k j 



A V 
MOM 



We get 



A^ = 0.34 ±0.05 GeV. (27) 



MS 

This is the main result of our computation. 

Our result can be directly compared with the one of ref. 



= 0.293 ± 0.018+°;°| GeV. (28) 

In terms of our result yields: 

a^(2.0 GeV) = 0.25 ± 0.02. (29) 

We do not attempt to estimate a^j for nj ^ on the basis of quenched 
data; a computation in full QCD with two degenerate flavours of sea quarks 
is in progress. 
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4.1 Starting from the Lattice 



One can apply renormalisation group considerations, analogous to the ones 
that establish the momentum dependence of the renormalised coupling, to the 
bare theory. On the lattice, as in any other regularisation scheme, renormal- 
isability implies that the bare coupling constant should be cut-off dependent 
and, for a small enough, g = go(a) should be a universal function of a. Thus 
in the scaling region one can define, up to an arbitrary integration constant, 
an a-independent Klatt parameter, in terms of the lattice spacing and of the 
bare lattice coupling as 

1 / ^0 W dg' \ 
Klatt = -exp(j — ). (30) 

In the asymptotic scaling region a — > and go(a) — ► 0, where the /3-function 
is perturbatively computable, one can fix the integration constant by defining 

Klatt as 

A^Ie^-^-i^) (31) 

independent of a in the asymptotic scaling limit. Note that, whereas we 
find asymptotic scaling for the renormalised coupling for momenta 1.8 < 



q < 2.3 GeV, see eq. (|2l|) and Figure 3, this does not necessarily imply that 
the asymptotic regime has already set in for the bare coupling used in the 
simulation. As a check, which involves lattice perturbation theory, we can 

1. take the non-perturbative determination of ^^q M -, 

2. compute the ratio Klatt j K^- M in lattice perturbation theory 

3. extract Klatt from 

a a Klatt , , 

Klatt = K~ M — — (32) 

MOM 

and compare the value obtained from eq. (|32|) with the one from eq. (|3lf) . 
The ratio 



A 



A 



1/1 1 



26 ° 9 2 MOM^ 2o 2 W 



+ 0{gl) (33) 
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has been calculated at one-loop by Hasenfratz and Hasenfratz |nj in the 
Feynman gauge A = 1, where A here is the gauge parameter. They found 

69.4. (34) 



Klatt 

We have checked their result by working in a general covariant gauge. The 
details of the LPTH calculation are given in the appendix. In the Landau 
gauge, we get 



KlATT 



54.6. (35) 



The number that we need to insert in (^) is not quite the above one, be- 
cause of the already mentioned extra term in the vertex renormalisation, not 
proportional to the tree level vertex. In our scheme the result is 



A mW A = °) 



83.2. (36) 



Klatt 

By inserting it into eq. fl32|), and using the result (|2|), one would get 

Klatt = 11.6 MeV, (37) 



Klatt = 4.5 MeV (38) 



to be compared with 



obtained from the hypothesis of asymptotic scaling, eq. fl3~T|), with g\ = 
1 and a -1 = 1.9 GeV. The comparison confirms the well known result 
that for values of f3 accessible to current simulations, Klatt still displays (3 
dependence 0. We stress again that in our case the matching procedure, 
described in the previous subsection, does not require knowledge of Klatt. 

Another way of seeing the failure of LPTH in this case is to compare the 
non-perturbative results for g 2 (/x) with what is obtained by inserting in the 
relation 

9 2 (t*) = Z\{^a)Zy 2 {^a)gl{a) (39) 

the values of the Z's obtained in LPTH and "boosted" lattice perturbation 
theory |J. This comparison is shown in Figure 5. 
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11 - 

10 - 
9 - 



Figure 5: g 2 (fi) vs. momentum as evaluated non-perturbatively on the 16 4 
lattice and from standard (solid line) and boosted (dashed line) LPTH. 



Finally, we can perform a cross-check of our perturbative calculations 
on the lattice and in the continuum as follows. It is known |18| that at the 
one-loop level 



A' 



MS 



Ai 



28.8. 



(40) 



This is a gauge-invariant quantity, since A^g^ and Klatt are both gauge 
independent. We can insert in the identity 



^-MS \(^MOM^'\ _ ^MS 



( 1V MS \ 

va — m/ 



MOM 



(A)- 



Ai 



(41) 



our continuum result for Aj^/A, 



and the lattice one for A ,7?^, ,j Alatt , 



MS/ AX MOM ""^ J -" ,UUJ -^ x ''MOM/ 

which we have evaluated for a generic covariant gauge. For any value of A 



we get indeed the value fl40|) found in the literature. 



5 Conclusions 



We have shown in the pure gauge theory that a non-perturbative determina- 
tion of the QCD running coupling can be obtained from first principles by a 
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lattice study of the triple-gluon vertex. We have some evidence that system- 
atic lattice effects are under control in our calculation. The main features of 
our method are that LPTH is not needed to match our results to MS and 
that the extension to the full theory does not present in principle any extra 
problem. Encouraged by our results in the quenched approximation, we are 
now repeating our calculation with dynamical quarks. 
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7 Appendix 



We start by giving the details of the (quenched) one-loop continuum per- 
turbative calculation needed to relate the scale parameter A in the MOM 
scheme to the one in MS. 



The three-gluon vertex in the MS scheme, calculated at the asymmetric 
point and at p 2 = fi 2 , can be written in a generic covariant gauge as (we 
suppress colour indices) 



r Ss a/37 GV/^°>-p/^)ip 2 =M 2 



+ % 



where A is the gauge parameter and 



CV(A) : 

The (Euclidean) tree-level vertex is 



1 + 



167P 



61 



167T 



-,C V {\) 



rCextra(A) 



ri-L af3~/(P> 0' P) 
PaPpP-f 



8 ai pp 



P 



(42) 



I8 + 2 A ' 
37 



6 



3A + -A' 



,(3) 

tree a/87 



(p, 0, -p) = -ig [5 a(3 p^ + 5 lP p a - 25 aj p /3 ) 



(43) 



(44) 



Note that despite the vanishing of one of the external momenta, the renor- 
malised vertex (021) is a finite quantity, as the infrared behaviour is completely 



controlled by the off-shell gluon | l?fl . 

Setting a = 7 and summing over a we get 



eigpp 



167T 2 



C V (X) + 



Cerrfra(A) 



(45) 



Whereas in the MS scheme only the pole part appears in the definition of 
renormalisation constants, nontrivial finite terms are included in momentum 
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subtraction schemes. In particular, in our MOM scheme the finite combina- 
tion (Cy + C ' ex tra /2), which appears in eq. fl45D , enters in the definition (0) 
of Zy. Hence, by computing 



Zl(e,^)Z v 2 (e,fi) gl(e) 



in the two schemes, one finds that the coefficient C- 
the two couplings according to (| 



MS, MOM' 

is given by the expression 



a 



MS, MOM 



where 



1 

47T 

C A {\) 



3C A (A) - 2 CV(A) + 



C extra (A) 



97 

36 



1a + Ia» 

2 4 



(46) 

which relates 
(47) 
(48) 



is the finite contribution coming from the gluon self-energy at one-loop [|T^ 
By setting A = 0, we obtain the result in the Landau gauge 



-1.856807669. 



MS, MOM 

Finally, the ratio (|25]) of the A parameters is obtained from 



MS 



A 



'MOM 



cxp 



n 

MS,MOM 

8tt b n 



0.35, 



(49) 



(50) 



where bn 



r-pj- in the quenched approximation. 



Turning now to our lattice calculations, we observe that the part of the 
one-loop three-gluon vertex on the lattice that is proportional to the tree-level 
vertex can be written as 



,(3) 



r (3) x 

1 tree A 

r (3) x 

1 trcc X 



1 + 9oA L 
(l + 9oA 



A/5 



+ 9oCi 



(51) 



where A-^ and A-j^ stand for the (momentum dependent) lattice and con- 
tinuum (MS scheme) contributions respectively. The quantity C*l, which 
relates the pure lattice result to the MS scheme, is momentum independent. 
Extra terms not proportional to the tree-level vertex are equal on the lat- 
tice and in the continuum, thus they do not contribute to C*l, which can be 
computed by evaluating the diagrams of Figure 6. 
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Figure 6: Diagrams for the one- loop three-gluon vertex on the lattice. Dashed 
and solid lines indicate ghosts and gluons respectively. 

The contribution of each one to Cx is 



a) 
b) 

c) 
d) 



f) 






N ( 1 r 1 j 1„ 

— t;L --J+ -Z 

24 V16vr 2 4 6 

N 

T 



16tt 2 4 16 



143 Z 1 



3 7 + 3 7 13 7 



1 

8iV 
AT 



13 



6l + 8 Z ° 



A)iV^o 



167T 2 4 



13 



16 



96tt ; 



; 8 V16tt 2 4 24 



+ 



(52) 



In these equations L, Z and J are [19, 20 



L = In /x 2 a 2 + 7euicr — In 4vr 
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, ^ t^4t4 = 0.1549333902311, 
I-* (2n) 4 q 2 

J = 0.0465621749414. (53) 

Finally /x is the mass scale introduced during the dimensional regularisation 
of the continuum MS scheme. These formulae are valid for the gauge group 
SU(N). Notice that both A L and Ajjg depend on A, A 2 and A 3 , while Cl 
depends only on A. 

Collecting all contributions in eq. fl52D , we obtain the result for Cl 



n 1 2 / 1 1 11 3 1 \ 

° L ~ m + 3 N {lQ^ L ~A J ~96 Z °~ 32 + 32^) + 

(l-X)-N(^—L--J + —Z ). (54) 
v U V16tt 2 4 24 V v ' 



Now, using the results of reference |l7j we get the result for Al in the 
Landau gauge 

A l = ?T7^ ln A 2 - 0.294728. (55) 

4 167T^ 
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